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ABSTRACT 



Context. Intriguing features in the angular distribution of the cosmic microwave background (CMB), such as the north-south asymme- 
try, were reported in the one- and three-year Wilkinson Microwave Anisotropy Probe (WMAP) data and should be studied in detail. 
We investigate some of these asymmetries in the CMB temperature angular distribution considering the ACDM model in the three, 
five and seven year WMAP data. 

Aims. We aim to analyze the four quadrants of the internal linear combination (ILC) CMB maps using three different Galactic cuts: 
the WMAP KQ85 mask, a \b\ < 10° Galactic cut, and the WMAP KQ85 mask + \b\ < 10° Galactic cut. 

Methods. We used the two-point angular correlation function (TPCF) in the WMAP maps for each of their quadrants. The same 
procedure was performed for 1000 Monte Carlo (MC) simulations that were produced using the WMAP team ACDM best-fit power 
spectrum. In addition, we changed the quadrupole and octopole amplitudes obtained from the ACDM model spectrum. We changed 
this to fit the quadrupole and octopole amplitudes to their observable values from the WMAP data. We repeated the analysis for the 
1000 simulations of this modified ACDM model, hereafter MACDM. 

Results. Our analysis showed asymmetries between the southeastern quadrant (SEQ) and the other quadrants (southwestern quad- 
rant (SWQ), northeastern quadrant (NEQ) and northwestern quadrant (NWQ)). Over all WMAP ILC maps, the probability for the 
occurrence of the SEQ-NEQ, SEQ-SWQ and SEQ-NWQ asymmetries varies from 0.1% (SEQ-NEQ) to 8.5% (SEQ-SWQ) using the 
KQ85 mask and the KQ85 mask + \b\ < 10° Galactic cut, respectively. We also calculated the probabilities for the MACDM using 
only the KQ85 mask and found no significant differences in the results. Moreover, the cold spot region located in the SEQ quadrant 
was covered with masks of 5,10 and 15 degrees radius and again the results remained unchanged. Furthermore, this analysis was re- 
peated for random regions in the SEQ quadrant with a 15-degree mask and the SEQ quadrant still remained asymmetric with respect 
to the other quadrants of the CMB map. 

Conclusions. We found an excess of power in the TPCF at scales > 100 degrees in the SEQ with respect to the other quadrants that 
is independent of the Galactic cut used. Moreover, we tested a possible relation between the Cold Spot and the SEQ excess of power 
and found no evidence for it. Finally, we could not find any specific region within the SEQ that might be considered responsible for 
the quadrant asymmetry. 

Key words. Cosmic Microwave Background - cosmology: observations - methods: data analysis - methods: statistical 



1. Introduction 

After the cosmic micro wave background (CMB) discovery by 
A. Penzias e R. Wilson dPenzias & WilsonL fl965l) . several ex- 
periments were developed to characterize this radiation, lead- 
ing to high-precision observations that raised the status of the 
CMB, which is since considered to be one of the main pillars of 
the ACDM model. Particularly, the CMB cosmological fluctua- 
tions first det ected by the Cosmic Background Explorer (COBE) 
satellite data dSmoot et al.LI 1992b set up a new era in cosmology 
studies and paved the way to what today is called precision cos- 
mology. However, detailed studies of the angular distribution of 
temperature fluctuations showed unexpected results if analyzed 
within the framework of the so-called cosmological concordance 
model (ACDM model). These peculiar features in the CMB an- 
gular distribution were found for the first time in the COBE data 
and attracted much interest since then. 

A quadrupole amplitude smaller than that expected ac- 
cording to the ACDM model was reported by the COBE team 



dSmootet all 1 19921) an d was c onfirmed by al l WMAP data 
releases dBennett et all 12003b iHinshaw et all I2007L 12009b 
iJarosik et all 1201 ll) . Other anomalies were found in the WMAP 
data that were not expected according to the ACDM model ei- 
ther, s uch as the alignment be t ween the quadrupole and octopole 
(e.g dBielewicz et all 12004 ISchwarz et all l2004t ICqpi et al " 



20041; Ide Oliveira-Costa et all 12004 iBielewicz et all 12005 



Land & MagueiioL \20oK ICopietali 
120061 iFrommert & EnsslinL 120101: 



20061 lAbramoetal 



iGruppuso & Burigam , 



2010 j)), the low qu ad rupole a nd octopole a mplitu de s (e.g. . 



dMukheriee & Wand. 120031: lAvata et all 1201 Ot ICavonL 
120101: Cruz et alT" 201oT)). the north-south asymmetry (e.g., 



dEriksen et all l2004t lHansen et all 120041: lEriksen et all 12004b : 
[ Hansen et all l2004bt [Donoghue et aTf l2005t iHoftuft et al. . 
20091: iPacietall 1201 Ot iPietrobon et all 1201 Ot IVielvaeta" . 
|2010|) ). the anomalous alignment of the CMB features towa rd 



the Ecliptic poles (e.g., dWiaux et all 120061 1 Vielva et alll2007l)) 
and the cold spot (e.g ., dVielva et all 12004 ICruz et all 120051 
[20071: IVielva.Pl 12010) ). 
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Recently, lAluri & Jainl d201 ll) anal yzed in detail the signa- 
ture of parity asymmetry first found by iKim & Naselskvl d2010h 
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in the WMAP best- fit temperat ure power spectrum, confirming 
this asymmetry on a 3-cr level. lAluri & Jainl (120111) also con- 
cluded that their result is not due to residual foregrounds or to 
foreground cleaning. In addition, the preferred direction of the 
parity asymmetry coincides with the CMB kinematic dipole, 
showing that it ma y somehow be related to the quadrupole- 
octopole alignment (Na selskv et all [201 ll) . 

On the other hand. lBennett et all (1201 ll) reviewed the anoma- 
lies reported in the CMB temperature fluctuations and claimed 
that they are not statistically significant and, for this reason, do 
not in disagree with the ACDM concordance model. 

Nevertheless, in this work we report an asymmetry in the 
WMAP temperature anisotropy data appearing in the two-point 
angular correlation function (TPCF) at scales above 100 degrees. 
In Section 2 we present our method to prepare the MC sky 
map simulations to confront them with real data and describe 
our TPCF method. In Section 3 we show our results. Finally, in 
Section 4 we present the discussion, followed by the conclusions 
in Section 5. 



2. Method 

The results presented in this paper were derived from the anal- 
ysis of three temperature ILC maps fr om the third, fifth an d 

" 2p07h . 



seventh year of WMAP data: WILC3 dHii 
WILC5 fcold et al.L[2008b and WILC7 (iJarosik et al.ll2011l) . We 



computed the TPCF for the WMAP data and for the MC simula- 
tions to obtain, evaluate and finally compare the results with the 
ACDM model. The HEALPix (hierarc hical equal area and iso- 
latitude pixelization) package (synfast) dGorski et al.Ll2005h was 
used to generate the MC simulations and analyze the maps. 

The simulated sky maps were generated with two different 
seed spectra. In the first run, we used the WMAP5 best-fit spec- 
trum to the ACDM model, available at LAMBDA^ to generate 
1000 simulations with N s id e = 256 (pixel diameter ~ 14'). In 
the second run, we modified the ACDM model spectrum, sub- 
stituting the amplitudes of the best-fit quadrupole and octopole 
by the values reported in the WMAP five-year data. This mod- 
ified spectrum was used to generate 1000 simulations with the 
same resolution as the first run (the MACDM model). 

The data and the simulations were analyzed by means of 
the TPCF. Since it is hard to analyze the CMB angular distribu- 
tion because of the foreground contamination, even using mul- 
tifrequency maps, we chose to use three different Galactic cuts. 
When the first CMB ma ps covering the whole sky we re released 
from the Relict satellite (IStrukov & Skulachevill984 . the safest 
way to deal with the foreground contamination was making par- 
allel cuts above and below the Galactic plane. More recently, the 
masks developed by the WMAP team remove in a fair way the 
known point sources and the Galactic signal, but they are not 
perfect since the Galactic foregrounds are still not completely 
known. [S moot et al] d 19921) also warned that the TPCF is highly 
affected when regions \b\ < 10° are included in the analyses. 
Taking into account these results, we used the WMAP KQ85 
mask, a Galactic cut \b\ < 10° and the WMAP KQ85 mask + 
\b\ < 10° Galactic cut in the present work. This last cut avoids 
the Galactic foreground and also the known point sources. All 
calculations were performed for degraded maps with N S id e = 64. 
We analyzed the NWQ with 9906, 10200 and 9495 pixels for 
each Galactic cut. For the SWQ, 10423, 10200 and 9795 pix- 
els were used, for the NEQ, 10250,10176 and 9417 pixels were 
used, and for the SEQ 9854, 101 76, 9350 pixels were used, when 



the Galactic foreground was removed with the WMAP KQ85 
mask, the Galactic cut \b\ < 10°, and the WMAP KQ85 mask + 
\b\ < 10° Galactic cut, respectively. We also counted the pixels 
left in each quadrant when the new WMAP mask KQ85y7 was 
used to remove the Galactic foregorund, there were 957 1 pixels 
in the NWQ, 10093 pixels in the SWQ, 9876 pixels in the NEQ, 
and 9217 pixels in the SEQ. 

The influence of the KQ85 mask asymmetry on the results 
was evaluated considering the three Galactic cuts. The same 
procedure was applied to the ACDM simulated sky maps. The 
MACDM realizations were analyzed using only the WMAP 
KQ85 mask. This modified spectrum enabled us to evaluate if 
the observed low quadrupole and octopole values can account 
for the results. We did not repeat the analysis using the KQ85y7 
mask for the MC simulations because the computational time for 
computing the TPCF is prohibitive. 

2. 1 . Two-point angular correlation function 

The TPCF function measures the angular correlation of temper- 
ature fluctuations distributed in the sky and is defined as 



c( r ) = <r(n p )r(n q )>. 



(1) 



7\iip) and 7\iiq) are the temperature fluctuations of the p and 
q pixels, respectively, and y is the angular distance between the 
two pixels. 

The pixels p and q are defined by the coordinates (0 p , <p p ) 
and (6 q , <f> q ), where 0° < (f> < 360° and -90° < 9 < 90°. It is now 
possible to obtain the equation for y: 



cos y - cos 6 p cos G q + sin 6 p sin 6 q cos((f> p - <p q ). 



(2) 



Finally, we define an rms-like quantity, cr, to compare the 
TPCF computed b oth for WMAP data and MC simulations 
dBernui et alll2006h : 



\ Nbins |} 



(3) 



The corresponds to the TPCF for each bin i. We used the 
number of bins Nbins = 90 to quantify our results. 



3. Results 

As mentioned before, we divided the CMB sky into quadrants 
and computed the TPCF for each quadrant, using the three dif- 
ferent Galactic cuts: the WMAP KQ85 mask, the Galactic cut 
\b\ < 10° and the WMAP KQ85 mask + \b\ < 10° Galactic cut. 
The results are shown in Figures Q] to [6] Asymmetries between 
the SEQ and the other quadrants can be easily noticed in the 
curves. On the other hand, the TPCF curves for the NWQ, the 
NEQ and the SWQ are mostly inside the gray-shadowed area, 
which corresponds to the MC lcr interval. 

Tables[l][2]and[3]show the probability for these asymmetries 
to occur for the ILC maps, and for the MC maps derived from 
the ACDM model. Both consider the above mentioned Galactic 
cuts. 

These tables also contain two sets of probabilities, PI and 
P2. The probability PI refers to the chance that exactly the same 
asymmetries found in the WMAP sky maps appear in the MC 
simulations: 



http://lambda.gsfc.nasa.gov 



- ctseqI o-NWQ(simul.) > o- SE Q/o- NW Q(data) 
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Fig. 1. TPCF curves computed for the WILC7 map using the 
WMAP KQ85 mask, smoothed for illustration purposes. The 
NWQ (top) and the NEQ (bottom) appear as solid red lines. 
The shadow part depicts the standard deviation intervals for 1000 
simulated maps produced with the ACDM spectrum. 

- cr SE Q/cr SW Q(simul.) > o- S EQ/o~swQ(data) 

- o- SEQ /cr N EQ(simul.) > o- SEQ /o- NE Q(data). 

The probability P2 extends the range of PI, comparing the 
chance that the ratio between the asymmetries found in the SEQ 
and in any of the three quadrants (from the MC simulations) 
to exceed that from the ratio between the SEQ and one given 
quadrant. The ** in the expressions below apply for any of the 
three quadrants (** = NW or S W or NE). 

- o- SEQ /o-„ Q (simul.) > o- SEQ /o- NWQ (data) 

- a- SEQ lcr„Q{simul.) > o- SEQ /o- SWQ (data) 

- o- SEQ /cr tt Q(simul.) > o- SEQ /o- NEQ (data). 

3. 1 . The cold spot 

Because the asymmetry was coincidentally found in th e same 
quadrant as the s o-called cold spot (see dCruz et al 1 l2007t 
IVielva et all 120101) '). we tested if they are possibly correlated. 
We masked the cold spot region which is centered on <p — 209° 
and = 141°, and calculated the TPCF for the SEQ using also 
the KQ85 mask. The cold spot mask radii were 5°, 10° and 15°. 

Finally, to look for a specific region in the SEQ that could 
account for the asymmetry, we scanned this quadrant, masking 



Fig. 2. Same as Figure [TJ but now with the SWQ at the top and 
the SEQ at the bottom. 



four different regions in addition to the cold spot region. The 
masks were centered on coordinates chosen randomly on (0 = 
315°, = 157°), (<p = 225°, = 113°), {<p = 270°, = 135°), 
(0 = 315°, = 113°), all of them with a radius of 15°. The 
TPCF was computed using a mask in each position at a time. 
The results are shown in Figure [7] (bottom). 

4. Discussion 

Our analyses indicate that, even though the excess of power in 
the SEQ is independent of the chosen Galactic cut, the largest 
asymmetry between the SEQ and the other quadrants depends on 
the Galactic cut and also on which ILC map is used. Considering 
the analysis performed only with the KQ85 mask, the biggest 
asymmetry was found between the SEQ and the NEQ for all 
three ILC maps, the largest corresponding to the WILC5 map. 
Moreover, when only the KQ85 mask was used, the asymme- 
tries between the SEQ and the other three quadrants for the other 
Galactic cuts are smaller. However, we found that the chance of 
having any asymmetry is quite low. Table 4 shows that most of 
the sigma ratios between the SEQ and the other three quadrants 
obtained from the TPCF computed from the data (see Eq. [3} are 
outside the standard deviation interval of the sigma ratios ob- 
tained from the TPCF computed for the MC simulations. 

Indeed, most of the sigma ratios between the SEQ and any 
other quadrant obtained from the TPCF computed for the simu- 
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Fig. 3. Same as FigureQ] but now using the WMAP KQ85 mask 
+ \b\ < 10° Galactic cut in the temperature WMAP7 data. From 
top to bottom, the curves refer to the NWQ and the NEQ, re- 
spectively. 



lated sky maps average one, as expected, because we generated 
realizations with randomly distributed temperature fluctuations 
in the sky. 

Finally, it is possible to see that the excess of power in the 
SEQ is independent of the chosen Galactic cut. However, the 
chosen Galactic cut influences the size of this excess in compar- 
ison to the other quadrants. Furthermore, using the asymmetric 
mask does not change our results drastically. We also notice this 
by analyzing the effect of the new mask provided for the WMAP 
sseven-year data release (KQ85y7) in TPCF There is no clear 
difference in the results of the TPCF for each quadrant using 
the WILC7 + KQ85 and WILC7 + KQ85y7 data sets, as can be 
seen from Figures[8]and|9] The values found for the sigma ratios 
between the SEQ and the other quadrants in this last case were 
oseq/o~nwq = 4.6, o- SE q/o-swq = 4.7 and o- seq /ct NE q = 6.8 
(see TableQ]for comparison). 

Furthermore, using the MACDM, which includes the 
WMAP best-fit values for the quadrupole and octopole, instead 
of the ACDM to generate the MC simulations does not notice- 
ably change neither PI nor P2, as can be seen in Table [5] Based 
upon these results, we are confident to state that the low values 
for the quadrupole and octopole in the data are not responsible 
for this quadrant asymmetry. 



Fig. 4. Same as Figure[3] From top to bottom, the curves refer to 
the SWQ and the SEQ, respectively. 



Finally, no evidence for a relationship between the SEQ 
asymmetry and the cold spot was found. The excess of power 
in the SEQ was most evident above the angular distance of 100 
degrees. By masking the cold spot, we were able to identify 
small differences in the behavior of the TPCF, all of them re- 
maining within the angular distance of 100 degrees. The same 
result was found for the other randomly chosen regions in the 
SEQ quadrant and no evidence was found that can relate the re- 
ported asymmetry to any of these chosen regions in the SEQ. 



5. Conclusion 

We found a significant asymmetry between the SEQ and the 
other quadrants by considering the temperature WMAP ILC 
maps. We calculated the probability of occurrence for this asym- 
metry using MC simulations, and 1 out of 1000 simulations for 
the ACDM model corresponded to the SEQ-NEQ asymmetry 
found in the WILC7 using the KQ85 mask. We also showed that 
different Galactic cuts do not influence the result in a significant 
way, leading us to believe that this effect is not caused by the 
asymmetric mask. Moreover, the use of KQ85y7 preserves the 
same asymmetries (SEQ-NEQ, SEQ-SWQ and SEQ-NWQ), as 
expected. Considering all Galactic cuts and maps used, the high- 
est probability of having an asymmetry is 8.5%, for the WILC5 
with the KQ85 mask + data clipping for \b\ < 10°. 
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Fig. 5. Same as Figures QJU but now using \b\ < 10° Galactic 
cut, again in the temperature WILC7 map. From top to bottom 
the curves refer to the NWQ and the NEQ, respectively. 



The possibility that the asymmetries described in this work 
were related to the reported lack of power in the quadrupole 
and octopole was tested. We constructed simulations based 
on a modified ACDM model, adjusting the amplitudes of the 
quadrupole and the octopole to their observational WMAP val- 
ues. No explicit relation between the quadrant asymmetries and 
the low amplitude of the first two non-zero multipoles was 
found. 

Furthermore, we found no evidence of a relationship be- 
tween the cold spot region an d the SEQ excess of power, as 
pointed out by iBernui who suggested that the Cold 

Spot is responsible for 60% of the Southern Hemisphere power. 
Masking this region and some other regions in the SEQ does not 
change the TPCF in a significant way, leading us to conclude 
that the asymmetries between the SEQ and the other quadrants 
are not related to any specific region in the SEQ. 

We conclude that the excess of power found in the SEQ is 
likely related to the north-south asymmetry, in which the South 
Hemisphere presents more power than the N orthern one (see 
dEriksen et ail 12004 lHansen et all [20041) and dPaci et al.Ll2010l 



Fig. 6. Same as Figure for the SWQ (top) and the SEQ (bot- 
tom). 

An explanation for these asymmetries is still missing. They 
could be primordial or caused by residual foregrounds or sys- 
tematic effects. The upcoming CMB data from the Planck satel- 
lite when analyzed with more accurate foreground removal tech- 
niques will enable us to study in more detail the CMB anomalies 
reported in the literature. 

Finally, in addition to the SEQ excess of power, we can no- 
tice a lack of correlation in the SWQ, NWQ and NEQ in the 
TPCF in scales between 20 and 100 degrees for all Galactic cuts 
in the present work (see Figures [T]|6). A lack of c orrelation in the 
TPCF was already reported bv lCopi et al.l d2007l) in scales above 
60 degrees for a full sky analysis using the KpO mask in the first 
and third year of WMAP data. 

We would like to thank Paolo Cabella for useful discus- 
sions. We also acknowledge the use of HEALPix packages, 
of the Legacy Archive for Microwave Background Data anal- 
ysis (LAMBDA). L. Santos thanks CAPES-Brazil for financial 
support. T. Villela acknowledges CNPq support through grant 
308113/2010-1. C. A. Wuensche acknowledges CNPq grant 
308202/2010-4. 
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planation for the north-south asymmetry. 
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The red dashed-dot curve refers to the Cold Spot region. 
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Fig. 9. Same as Figure |8] for the SWQ (top) and the SEQ (bot- 
tom). 
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Table 1. Calculated probabilities of finding the same asymmetries as in WMAP data in the MC simulations using the WMAP KQ85 
mask and considering the ACDM model. 
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(a) Probability of finding, in this case, the asymmetry between the SEQ and the NWQ quadrants in the simulations. {b> Probability of finding the 
asymmetry between the SEQ quadrant and any other quadrant in the simulations. 

Table 2. Calculated probabilities of finding the same asymmetries as in WMAP data in the MC simulations using the WMAP KQ85 
mask + \b\ < 10° Galactic cut and considering the ACDM model. 



Map 


c seq/cnwq 


PI 


o~seqIcswq 


PI 


o~seq/o~neq 


PI 


WILC7 


2.7 


4.9% 


2.3 


7.4% 


2.5 


6.0% 


WILC5 


2.3 


7.4% 


2.2 


8.5% 


3.2 


2.7% 


WILC3 


2.4 


6.5% 


2.4 


6.5% 


3.0 


3.1% 


Map 


°" seq/cnwq 


P2 


cr SEQ/ cr SWQ 


P2 


o~seq/o~neq 


P2 


WILC7 


2.7 


10.1% 


2.3 


16.0% 


2.5 


12.8% 


WILC5 


2.3 


16.0% 


2.2 


18.3% 


3.2 


5.4% 


WILC3 


2.4 


14.0% 


2.4 


14.0% 


3.0 


6.5% 



Table 3. Calculated probabilities of finding the same asymmetries as in WMAP data in the MC simulations using the \b\ < 10° 
Galactic cut and considering the ACDM model. 



Map 


°" seq/cnwq 


PI 


cr SEQl cr SWQ 


PI 


Cseq/cneq 


PI 


WILC7 


3.7 


1.6% 


4.0 


1.2% 


2.8 


4.8% 


WILC5 


3.6 


1.8% 


3.6 


1.8% 


3.4 


2.5% 


WILC3 


3.4 


2.5% 


3.6 


1.8% 


3.1 


3.2% 


Map 


c seq/o~nwq 


P2 


Cseq/ctswq 


P2 


o~seq/o~neq 


P2 


WILC7 


3.7 


3.2% 


4.0 


2.8% 


2.8 


8.8% 


WILC5 


3.6 


3.4% 


3.6 


3.4% 


3.4 


4.1% 


WILC3 


3.4 


4.1% 


3.6 


3.4% 


3.1 


5.8% 



Table 4. Sigma ratios for the simulated CMB maps considering the ACDM model and the three Galactic cuts. 



Galactic cut 


Cseq/o'nwq 


CseqIcswq 


Cseq/cneq 


KQ85 


1.0 +LI 


1.1 +L2 


0.9 +10 


KQ85 + \b\ < 10° 


1.0+H 


i.o + H 


l.o+E 


\b\ < 10° 


1 +< ^ 


1 +< ^ 


1 o + ^ 



Table 5. Calculated probabilities of finding the same asymmetries as in WMAP data in the MC simulations using the WMAP KQ85 
mask and considering the MACDM model. 



Map 


c seq/o~nwq 


PI 


c seq/o~swq 


PI 


o~seq/o~neq 


PI 


WILC7 


4.6 


0.7% 


4.6 


0.7% 


6.7 


0.2% 


WILC5 


4.2 


1.3% 


4.1 


1.3% 


7.0 


0.2% 


WILC3 


3.9 


1.5% 


4.1 


1.3% 


6.8 


0.2% 


Map 


°" seq/o~nwq 


P2 


Cseq/cs wq 


P2 


Cseq/cneq 


P2 


WILC7 


4.6 


1.7% 


4.6 


1.7% 


6.7 


0.4% 


WILC5 


4.2 


2.6% 


4.1 


2.9% 


7.0 


0.3% 


WILC3 


3.9 


3.9% 


4.1 


2.9% 


6.8 


0.4% 
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